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Abstract. The aim of this paper is to describe efficient algorithms for computing Maass 
waveforms on subgroups of the modular group PSL{2,Z) with general multiplier systems and 
real weight. A selection of numerical results obtained with these algorithms is also presented. 
Certain operators acting on the spaces of interest are also discussed. The specific phenom- 
ena that were investigated include the Shimura correspondence for Maass waveforms and the 
behavior of the weight-fc Laplace spectra for the modular surface as the weight approaches 0. 



1. Introduction and notation 

The purpose of the present paper is to present computational methods and experimental 
results for Maass waveforms with general real weight and general multiplier systems. 

The classical theory of holomorphic automorphic forms was developed in the setting of (even) 
integer weight. This was motivated both from a geometrical point of view and by number 
theoretical applications (e.g. the study of modular forms related to the modular invariant j and 
the discriminant A). The problem of finding the number of representations of an integer as a 
sum of a fixed number of squares (cf. e.g. [l5] and [Hj) was successfully treated by using the 
theory of half- integral weight forms (e.g. 0-series). This motivated Petersson to develop a theory 
of automorphic forms and multiplier systems of arbitrary real weight [52] and later also complex 
weight [H I-IV]. 

Many applications of modular forms use Hecke operators as a principal tool and Wohlfahrt 
developed a theory of Hecke-like operators in arbitrary real weight [76] (cf. also [70]) . 

Now recall the definition of a Maass waveform as a real-analytic square-integrable eigenfunc- 
tion of the Laplacian on a Riemann surface of finite volume with constant negative curvature 
— 1 . These were introduced by Maass for zero weight in [40] . In [41] the theory of waveforms for 
general weights was developed by using the so called lowering and raising operators, which send 
a waveform of a given weight to one with a smaller or larger weight. 

In a sHghtly different setting Selberg [58l pp. 82-83] observed that considering the invariant 
differential operators on the space 5{ x S'^ with representation x and separating variables lead 
to real analytic eigenfunctions of the form / (z, 0) = F (z) e~'^'"^ where F (z) is a holomorphic 
modular form of integer weight k and character x- 

For an overview of the spectral theory of real analytic modular forms with arbitrary real 
weights and multiplier systems see e.g. Maass [12] and Roelcke [5^. More recent work can be 
found in e.g. [iOl [HOl [51] . [221123] and [SlIOlfTO]. 

It is also worth mentioning that the recent interest in Maass waveforms as representative 
objects for studying quantum chaos also applies to real weights. If a weight zero waveform 
corresponds to a quantum mechanical particle moving freely on a Riemann surface then a weight 
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k waveform represents a similar particle moving in a constant magnetic field with field strength 
proportional to k. 

1.1. Algorithms 

Previously published algorithms for computing Maass waveforms on cofinite Fuchsian groups 
have been restricted to groups with one cusp, e.g. the full modular group or Hecke triangle 
groups (see e.g. |63 [ [75 j [25 j [26 t [27 t l46j). By adjoining certain elements it is also possible to bring 
certain Hecke congruence subgroups to the one-cusp case, cf. e.g. [TB I fTT t l4]. 

In addition, with the exception of the (somewhat crude) computations in [36], only trivial 
multiplier system and zero (or even integer for the holomorphic case in [26]) weight has been 
considered. 

The most stable of the algorithms cited above is the one based on "implicit automorphy" 
by Hejhal (as detailed in e.g. [27]) which admits generalizations first of all to remarkably large 
spectral parameter (cf. [69]) and a further advantage is that it does not depend on any underlying 
arithmetical properties (i.e. Hecke operators). 

In [65 this algorithm was generaHzed to groups with several cusps, e.g. Hecke congruence 
subgroups ro(A^) with non-trivial Dirichlet characters and in and [6H ch. 3] we also considered 
general subgroups of the modular group. Recently, in [6], this algorithm in combination with 
other theoretical methods was used to show the existence of certain Maass waveforms close to 
the tentative waveforms produced by the algorithm. 

The aim of the present paper (which is based on [BU ch. 2]) is to demonstrate how to extend 
the algorithm to general multiplier systems and arbitrary real weights. 

The first section contains a brief review of the basic theory of multiplier systems and then we 
will introduce the notion of Maass waveforms in this context. We will also provide some details 
on the different operators that act on the space of Maass waveforms with non-trivial multipHer 
system. While being of interest in themselves, these operators can also be used in combination 
with other tests of reliability and accuracy of the algorithm. 

In section [7] we will give the specifics of how the algorithm is modified and in the last section 
we present a selection of the results which has been obtained with the described method. 

1.2. Summary of notation 

We will use the notation e{x) = e^'^*^ and for a complex number z we always use the principal 
branch of the argument, — tt < Argz < tt. 

Let Oi = {z = X + iy \ y > 0} he the upper half-plane equipped with the hyperbolic line- and 

area-elements ds^ = -^^^ and d/x = respectively. The boundary of J{ is = M U {oo}. 



The isometry group of J{ is identified with PGL(2,M) = GL(2, M)/ {±Id}, where Gi(2,R) is 
the group of invertible two-by-two matrices with real elements and Id = ( g i ) • For j — ( " ^ ) e 
GL(2,M) and 2 € IK we define an action by 



72: = < 



az+b 
cz+d ' 



if ad — be > 0, 



if ad-bc<0. 



The subgroup of orientation-preserving isometrics of "K is given by PSL{2, R) ~ SL{2, M) / {±Id} 
where SL{2,R) is the subgroup of GL{2,R) consisting of matrices with determinant 1. For any 
subgroup r C PGL(2,R), we use T to denote the inverse image of F in GL(2,R). Note that 
this forces —Id S F. We are mainly interested in Fuchsian groups, i.e. discrete subgroups of 
PSL{2,M.). Of particular interest is the subgroup consisting of matrices with integer entries, the 
modular group, PSL{2,Z). We are also interested in the so-called Hecke congruence subgroups, 
Fo(iV) = {(c d) G PSL{2,Z) |c = mod N} , defined for any positive integer N (note that 
Fo(l) =P^i(2,Z)). 
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We say that an element 7 of PSL{2, M) is elliptic, parabolic or hyperbolic if the absolute value 
of the trace of the associated matrix is smaller than, equal to or greater than 2 respectively, or 
equivalently, if 7 has one fixed point in IK, one (double) fixed point in 9IK or two (different) fixed 
points in Fixed points of paraboHc elements are called cusps. 

If r C PSL{2,M.) is a finitely generated Fuchsian group we identify the set of F-orbits with 
a connected subset of Jf, T =r\!K, a fundamental domain of F. If F has a set pi, . . . ,Pk of 
inequivalent cusps then the set T will meet dOi at k inequivalent points which we will also 
denote by pi , . . . , and we usually abuse the notation and call these points the cusps of F (or 
of T). By conjugation we may always assume pi ~ ioo. Corresponding to a cusp pj of F we use 
Sj to denote a parabolic generator of Tp. - the subgroup of F which fixes pj . We also choose a 
cusp normalizer, aj £ PSL{2,1S.), with the property that aj (00) = pj and ajSaJ^ = Sj, where 

= ( J J ) is the parabolic generator of PSL(2, Z) (the other generator being T ~ ( ? V ) • '^^^ 
map (Jj is uniquely determined up to a translation S'^ . H T meet at the points Qi, I < j < kq 
we fix a set of maps Ui ET such that Uiqi — pj where pj is the unique cusp equivalent to qi. 

2. Multiplier systems 

2.1. Introduction 

We will give a brief introduction to multiplier systems, for more extensive treatments see |22[ 
pp. 331-338], [53], or [55l pp. 70-87]. 

Let F be a Fuchsian group and m an even integer. Classically, a function ip, meromorphic on 
the upper half-plane IK, which satisfies 



(2.1) (p{Az)^eA{z;m)p{z)^{cz + d)"'p{z), VA^l " " ) e F, 



b 



c 



d 



is called an automorphic form of weight m for F. The function 

a b 



QA{z;m) = {cz + d)"\ A=\ ^ ^ )eF 



c 



is said to be an automorphy factor on F. The classical theory of automorphic forms is well-known; 
for instance, if m = 2, then the automorphic forms can be identified with the meromorphic 
differential forms of degree 1 on the orbifold (classical Riemann surface) F\IK. 

We observe that, for even m, the number {cz + c?)™ is uniquely defined and the automorphy 
factor 8^1(2:; m) in l|2.ip clearly satisfies 

(*) Qa{Bz; ■m)QB{z; m) = 6ab(z; m). 

To generalize these notions to arbitrary real m, there needs to be a choice of branch of the 
argument, and to make certain everything is well-defined, we have to introduce the notion of a 
multiplier system. 

Definition 2.1. For any real number m define 
To adapt the relation (*), we also write 

(2.2) (7m{A,B) = jA{Bz;m)jBiz;m)jAB{z;m)~^. 

It is clear that for integer m, cr,„(A, i?) = 1, but it can also be shown (cf. |53} §2, pp. 42-50]) 
that the only values which am can take are 1 and e*^'^"". 

Definition 2.2. v : T ^ ~ {z \ \z\ ~ 1} is said to be a multiplier system of weight m on F if 
• v{^I) = e-"™, and 
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• v{AB) = crrn{A,B)v{A)viB), yA,Be T. 

Observe that v can be regarded equally well as a multiplier system of any weight to' = 
m mod 2. The question of whether there exist multiplier systems of a given weight and on a 
given group is most easily answered by the following proposition (cf. [HI Prop. 2.1, p. 333]). 

Proposition 2.3. Given v :T ^ and m e M. The following are equivalent: 

• v{T) is a multiplier system of weight m on T. 

• There exists a function ^ on 3{ which is either C°° or meromorphic such that 

^(Az) = v{AMz)icz + dy\ yA^l^l ^ ^ e r. 

A (f a,s above is said to be an automorphic form of weight to and multipHer system z; on F. 

Definition 2.4. Given a multiplier system w on F and an element a G G'i(2,R) we define a 
multiplier system, w", on the group a~^Ta by 

am [a, A) 

That this indeed gives a multiplier system on a~^Ta is shown in [42l p. 138]. 

Using Prop. 12. 3| we will construct the two most widely used multiplier systems in the following 
sections. Compare: |22} pp. 334-337]. 
2.2. The rj multiplier system 

2.2.1. The ij-function. The Dedekind 77— function is a holomorphic function on "K, defined by 



00 

V{z) = e(^^)l[il-e{nz)). 



ri=l 

It is clear from the definition that 77(2) 7^ for z € 'K and that, for each € R, rf'^ can be defined 
as a holomorphic function on IK. (Cf. |55[ p. 205].) Note that ri{zY^ is the famous Discriminant 
function A(z). Cf. |55} pp. 196-197]. It is clear that ry (^ + 1) = e (^) and as for A [z) it is also 
possible to express rj{z) as a lacunary Fourier series (cf. [IHl p. 18]). 

2.2.2. The multiplier system. It can be proved (cf. Thm. 3.1 and Thm. 3.4 [H p. 48 and p. 52]) 
that 77 satisfies the following functional equations 

77 [ — ) = (—iz) 2 ry (z) , and in general. 



(2.3) 77(Az) = «,(A)(cz + d)H(z),VA= ^ ^)e5L(2,Z). 

This functional equation expresses the fact that 77 is an 5*^(2, Z)— automorphic form of weight \ 
and multiplier system given by Vj^. Accordingly the function rf'^ is an SL{2, Z)-automorphic form 
of weight k and multiplier system given by v'^'^ , and we can use rj'^'^ in the context of Proposition 
l2.3l to assure the existence of the multiplier system, Vjj^k = v^'^, of weight k on SL{2, Z) (and any 
of its subgroups, e.g. ro(A^)). We have the following explicit formula for v = v'^^: 

,^ 1 . (( a b\\ ( f|, a = d=l,c=0, 

^'■'^ ^^"^HU jj=| fc(^-.(rf,c)), 00, 

and for c < we use that v {—A) = e^'^^^v [A) (cf. Def. 12. 2p . Here s{d,c) is the Dedekind sum, 

n—1 ^ ^ ^ 
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where ((a;)) is the saw-tooth function 

x — [xj — i, if x ^ Z, 



((^)) 



if a; S Z, 



and [xj is the greatest integer less than or equal to x. Note that if x is not an integer, then 
[— a;J = — [xj — 1 so ((— x)) = — ((x)) , and hence s(— d, c) = —s{d, c) if gcd{d, c) = 1. 

Remark 2.5. It is also possible to express the eta multiplier explicitly without Dedekind sums 
but using extended quadratic residue symbols instead. We have the following formulas from 
Knopp [Ml P- 51] or van Lint [HI Thm. 3]: 



(2.5) 



a 6 \\ f (f) [(a + <i)c- M(c2 - 1) +3d- 3- 3cd]) , c> 0, even, 
c jj " 1 (l)e(^ [(a + d)c-M(c2-l) -3c]), c> 0, odd. 



(Note that the symbols (4)^ and (f)* of [3T|[71] agree with our symbols in these two cases.) 

Remark 2.6. It is known that, for each /c € K, there exist exactly 6 different multiplier systems 
of weight k on PSL{2,1) (cf. [55l §3.4, pp. 83, 206 ] or O Thm. 19, p. 132]). We will denote 
these by v^J^\ = v^l-'^^^\ where r € {0,2,4,6,8,10}. Compare [55l eq. (6.4.7)]; one knows, of 
course, that v^'^ ~ 1. 

When dealing with the modular group and weight k, it is sufficient to consider only the 
multiplier system v^l, = w^^fc = v^'^ (for reasons to be discussed later in Section [4.4. ip . 

2.3. The 9 multiplier system 

On any subgroup of PS'L(2,Z), we can always use the 77-multiplier system, but in general, 
on subgroups of P5'L(2,Z), there are also other multiplier systems available. In particular, on 
Fq (4) , there is a multiplier system of weight ^ which is interesting from an arithmetical point of 
view. 

It is well-known (cf. [60] or [28]) that the Jacobi theta function 

00 

0(^)=^e(n2z), 
—00 

is automorphic on ro(4) with weight fc = 5 and can be used to define a multiplier system on 
To (4). Using the Poisson summation formula one can prove (cf. [H pp. 72-75] or [28l pp. 
167-168]) that the theta function satisfies: 

(2.6) 

and one can also prove the general formula (cf. |28[ Thm. 10.10, p. 177] or [6OI p. 447]): 



(2.7) e{Az) = «,(A)(cz + d)t0(z),^= dj^ro(4). 

The multiplier vg (A) can be expressed explicitly as 

ve (A) = (^) , 

where = 1 if = 1 mod 4 and = * if ^ = — 1 mod 4, and (^) denotes the extended quadratic 
residue symbol defined as the traditional Jacobi symbol if < d = 1 mod 2 and extended by 



b 
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and 

'0\ fl ifd = ±l, 



J [0 otherwise. 

For the sake of completeness we also use the traditional Kronecker extension, i.e. we define 

One can verify that our symbol (7) satisfies reciprocity relations similar to the usual ones: 
Proposition 2.7. Suppose that c,d ^2, are odd and c 7^ 0. Then we have: 



-1 
~d 

c 



|_ (I) , d,andc<{). 



Remark 2.8. Relations l|2.6p and l|2.7p can also be proved using the corresponding relations l|2.3p 
for 77 and the following relation between the rj and the 9 functions (cf. [28l p. 177] or [311 Thm. 

12, p. 46]): e{z) = 

3. Maass waveforms 
The slash-operator — fi^z) can be extended to an operator of weight k as: 

f\lk,A]{z) = f{Az)jA{z;k)~'^, 

and the natural analog of the Laplace-Beltrami operator. A, which is invariant under this action 
is the weight-fc Laplacian: 

Afe = A - lyk— = y — ^ + — ^ - lyk- 



dx ' \dx^ dy^ J ' dx 

K r is a Fuchsian group we define the space M(r, d, k, A) consisting of Maass waveforms on F, 
of weight fc, multiplier system v and eigenvalue A, as the space of functions which satisfy the 
following conditions: 

1) /|[^,,](z) = KA)/(z),VA= ^ J ) er, 

2) Afc/ + A/ = 0, and 

3) \f\''d^l < cx). 

Observe that condition 1) is equivalent to 

1') fiAz) ^ viAluiz; k)f(z), VA = f |j ) e F. 



For purposes of the computational work to he described in this paper, we shall he content to 
restrict ourselves to cases where A > ^. (Cf. also here para. 4 of sect. \8.1\ helow.) 

Instead of the Bessel equation in the case of weight 0, condition 2) above gives us the Whittaker 
equation, and using the method of separation of variables gives us Whittaker functions instead 
of the K-Bessel functions at weight (for complete details see |22[ Chap. 9]). Since f{x + iy) is 
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no longer periodic in x, but instead satisfies f{z + 1) = v{S)f{z) = e{a)f{z), with a G [0, 1), the 
Fourier series of / can ([22l pp. 26, 348, 420(19)]) be written as 

oo , > 

(3-1) f{z) = ^f=^W^ssn(«+a)|,zfl(47r|n + «|2/)e((n + «)x), 

where Wi^fj,{x) is the Whittaker function in standard notation (cf. |15} vol. I, p. 264]) and R is 
the usual spectral parameter, A = | + R^. One notes here that Wo^iji{x) = Tr~ixi Km (|) . For 
k = 0, the expansion above thus reduces to usual Fourier expansion with 2yi Kiii{2Tr\n + a\y) as 
in [221 P- 26, prop. 4.12]. 

If we have more than one cusp we define functions fj related to / at each cusp, pj , of F by using 
the cusp normalizing maps aj from section [Ol and setting fj{z) = f\[aj.k] (z) = jaj {z; k)~^ f{ajz). 
It is easy to see that 

f,{z + l) = v{Sj)fj{z)^eia,)f,{z), 
with aj G [0, 1) (cf. |28[ p. 41]). Thus the Fourier series of / at the cusp j can be written as 



oo , s 



(3.2) fj{z) = 2^ / | W^sg»(n+a,)4,»fl(47r|n + aj-|y)e((n + aj». 



71 — — oo 



As in the case of weight and Dirichlet character, we say that the cusp number j is open or 
singular if aj = and closed if Uj ^ 0. If all cusps of F are singular for the multiplier system v 
we say that w is a singular multiplier system for F. 

Remark 3.1. Observe that, for the eta-multipHer on PSL{2, Z) and weight k, we have a = ai = 

k_ 
12 ■ 

3.1. Decomposition of the discrete spectrum 

It is known (see for example [TT] or [22]) that closed cusps (i.e. v{Sj) ^ 1) do not contribute 
to the continuous spectrum, and if all cusps are closed there is only the discrete part of the 
spectrum left and this is spanned by the Maass waveforms. We also know (see [22l P- 385]) that 
on the modular group with weight k the smallest eigenvalue is 

, _\Jif,_\Ji 

or larger. In the case of PSL{2, Z) and A: > 0, F{z) — y^ri{z)'^'^ has eigenvalue equal to Xmin- 

In this paper, any eigenvalues A G [Xmin, \] will be regarded as exceptional. The non- 
exceptional eigenvalues thus satisfy | < Aq < Ai < • ■ ■ < A„ — > oo. One can obtain lower bounds 
for the eigenvalue Aq (see [HI p. 183]), but in light of the numerical experiments in Section [ST 
they are not very effective (cf. Figure [21). 



4. Operators 

4.1. Conjugation and reflection 

Let J and K denote the refiection, Jz — — z and conjugation, Kz = z. Then J and K act as 
involutions on the space of Maass waveforms via the operations 



Kf = fiK{z)^f{z), 
Jf = f\j{z) = fi-z). 
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It follows from the definition of the action of GL{2,M.) on IK (cf. page[2|) that we can use the 
matrix J = (o -i) in Gi(2,R) to represent the operator J. For A= ^) we define 



A* = JAJ 



-1 



a —b 
— c d 



and then A** = A, and A{z)^j^ ~ —A*(z\ meaning that also —A{z) = A*{—z). 

Remark 4.1. It is easy to verify that if / € M{To{N),v, k, X), then Kf G M(ro(iV), tJ, -fc, A) 
and Jf G M(ro(A^), w*, — fc, A), where v* is the multipHer system determined by 

Of particular interest is the involution obtained by combining J and K, i.e. 

KJf{z) = f^jK{z)^J{^. 

It is easily seen that if / has Fourier coefficients Cj{n), then (by [HI p. 265(8)]) f\jK has Fourier 
coefficients Cj{n), and we thus would like to have / and f^jx belonging to the same space (i.e. 
transform according to the same multiplier system), since then we can assume that the Fourier 
coefficients are real. 

It is clear that if / G M{To{N),v, k, X), then KJf G M{To{N),v*,k, X) so we are left to see 
whether v* = v or not. By using the explicit formulas one can verify that indeed v* = v for both 
ve and (see |64[ pp. 66-68] for details) and we arrive at the following proposition. 

Proposition 4.2. If v is either the rj- or the 6-multiplier system (in the latter case '^\N) then 
a basis {pi, . . . , (7„i} o/ M(ro(A^), w, fc, A) can be chosen so that each gj can be expanded in a 
Fourier series at oo with real coefficients. 

Proof. We have seen that for both the theta and the eta multipHer systems the product KJ 
is a conjugate-linear involution of the space M{ro{N),v,k,X), and hence we can assume that 
any / G M.{To{N),v, k, X) is an eigenfunction of KJ with eigenvalue e, where |e| = 1. Note 
that if f{z) has a Fourier series expansion as above with Fourier coefficients c(n) then f\x,j has 
Fourier coefficients c(n) and hence c{n) ~ ec{n). Finally we observe that if e = e'^we can look 
at the function g = e'^/ which then satisfies KJg = e~^^KJf = e~*2e*^/ = g. After proper 
normalization it is thus no restriction to assume that the eigenvalue of KJ is e = 1, and that the 
Fourier coefficients are real. □ 

Remark 4.3. For the sake of completeness it should be remarked that in general one can not 
simultaneously take Fourier coefficients at cusps other than co to be real (cf. the next subsection 
where we introduce the map wat, which is a cusp normalizing map for the cusp at and which 
has eigenvalues itz"'^). 



4.2. The involution rjv 

As in the case of zero weight (e.g. [65]) we define ujnz = or equivalently ujn ~ ^ ^ 

We know that ujn is an involution of ro(A^), i.e. ro(iV) = ll;nTo{N)uj:^^ , but the question is how 
it relates to the weight and multipHer system. 

If / G M(To{N),v,k,X) it is easy to see that f\[k,uj,^] e M(ro(A^), w"", fc, A), and it is also 
easy to verify that t;"™(T) — v [lunTuj]^^) . We also have 

/|[fc,t^«]|[A;,c^„](2) = j„„(z;fc)~V('^Af2)|jj^^^^] = j^^{z;k)-'^j^^{ujNz;k)-'^f (ujjfz) 

_ ^~ikArg(VNz)^-ikArg{-l/\/Nz) j ^-^-j _ ^-i-rrk J^^)^ 
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and hence if we define the operator rjv by 
we have that 

= Id. 

To show that rjy is a linear involution, we also have to verify that u"" — v. This is easily done 
in the two cases iV = 1 together with v = Vri and iV = 4 together with v = vg. And we conclude: 

Proposition 4.4. If N = 1 and v = w,,, or N = A and v ~ ve the operator 

rjv :M(ro(iV),«,fc,A) ^M(ro(iV),f,fc,A), 

defined by 

is a linear involution, i.e. tn (af) ~ ar^f for all a ^ C and r'^f = f. Hence it has eigenvalues 
±1. 

Remark 4.5. Note that in the case of To (4) (i.e. iV = 4) if r^fiz) = ±f{z), then 

which means that the Fourier coefficients at the cusp at are proportional to the coefficients at 
zoo: 

(4.1) C2(n) = ±e-*^fcci(ri) = ±i-^ci{n). 
4.3. The operator L 

Definition 4.6. For iV = 4 and the 6'-multiplier system at weight ^, following Kohnen [32] or 
Katok-Sarnak [29j we define the operator L acting on M(ro(4), vg,^,R) by 

where T'"\ is the following (exceptional) Hecke operator 

^> 2 

(4.2) = ^ E / 

j mod 4 

Cf. Subsection l4.5.1l It can be shown that the operator L is self-adjoint, commutes with Ai and 
all Hecke operators T^^ ii P 7^ 2 (cf. I|4.9p ) and have the eigenvalues 1 and —5. Cf. [32] and [48j . 

Suppose that f{z) has Fourier expansions a la p.2p at the cusps pi — 00, P2 = and Pi, = — \ 
with Fourier coefficients in) respectively. By calculations similar to [32] or [32l [S] it can be 
shown that Lf has Fourier coefficients 6(n) (with respect to 00) given by 

{ai(n) + (l + i)a2(f), 71 = mod 4, 

al(n) + ^/2a3(^)(-l)'^, n = 1 mod 4, 
— ai(n), n= 2, 3 mod 4. 

Let C M(ro(4), Vg, \, A) denote the subspace introduced by Kohnen in [32], i.e. y"*" consists 
of ah / with Fourier coefficients ai(n) = for n = 2, 3 mod 4. By using l|4.3p it is easy to verify 
that is exactly the eigenspace of L corresponding to the eigenvalue 1. Unfortunately the 
eigenspace V~ corresponding to the eigenvalue — ^ is not as simple to characterize. However, 
l|4.3p can be used to identify V~ by means of certain relations between coefficients at the cusps 
at 00, and -~\. E.g. if ai(l) = 1 then 03(0) = -\/2 and if ai(l) = then 03(0) = 0. 




10 



FREDRIK STROMBERG 



Clearly does not in general commute with L but in case Lf = A/, if = A4/ and 

tnJ = e/ (with e = ±1) then ai (4) = e^/2X = — "\/2e. This should be compared with the 
results on newforms at weight zero, e.g. [U p. 147] and [HI p. 31 ]. 

4.4. Maass operators 

So far, the operators we have seen act on spaces of Maass waveforms of a given weight and 
multiplier system. 

We will show that we may Hmit the range of weights k to investigate to k £ [0,6]. For this 
we will use the Maass lowering and raising operators, i?^, which raise or lower the weight of a 
Maass waveform by units of 2. They are defined by 

rn± . d , d k 

and using the relation between the Whittaker function and the confiuent hypergeometric function 
together with the transformation formulas [iSl p. 258, (10)] (see also: [lU p. 183 (middle)] and 
[i3l p. 302 fines -3 and -1]) we see that for A: > (here we set Y = 47r|7i + a\y and tIq = n + a) 

(4.4) S+ (W^i ,^(r)e(n,.T)) ^ - Wt+2 ^^j^ {Y) e{n^x), n„ > 0, 

(4.5) E- (w. ^^jj{Y)e{n^x)'^ = - + i + ^^,^R e{n^^),n^ > 0, 

(4.6) E+ (w_^^,j,{Y)eM)= (^M^±^ + \ + ^') ^-^^,^B, (^) ^K^)' < 0' 

(4.7) E- {w_K^,RiY)e{n^x))=W_^^^j^ {Y) e(n„x), 71, < 0. 

To verify that they respect the weight fc-action it is easiest to proceed straight forward but to 
use the following form of the operators 

, -^ d k 

rn- , d k 

and write e-»^'-9(=^+'*) = ' . Cf. 021 p. 178]. 

One can then show that E^ maps M(ro(A^), v, k, A) into M(ro(A''), v,k±2, A), and that the 
composition 

E^^^E^ :M(ro(iV),t;,A:,A) K^M(ro(iV),«,fc,A) 

is just multiplication by a constant, which is readily seen to be nonzero anytime A > |. Hence 
E^ acts bijectively on the spaces corresponding to non-exceptional eigenvalues, i.e. they are 
always bijections for A > i. 

4.4.1. Maass operators and the symmetry about k = 6 \4-41\ First of all, observe that E^ only 
change the weight and not the multiplier system v, but in view of the remark after Def. 12.21 it 

.,2k 



is clear that 1; = is also a multiplier system of weight k + r for any r € 2Z, and with the 

t 

E+ : M(ro(7V), v^;l,k, A) M(ro(A^), 47+'^, A: + 2, A), 



notation v^l. = v'^^'^'^^'' it is clear that 



and 



E^ : M(ro(iV), v^;l, k, A) ^ M(ro(W), k - 2, A). 
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Our main purpose is to investigate the eigenvalues of Maass waveforms on the modular group 
when the weight and multiplier system are varied. That is, we would like to investigate the space 
M{ro{l),vl^l, k, A) for all A: e R and r G {0, 2, 4, 6, 8, 10}. Suppose that A > i so the the lowering 

and raising operators act bijectively, then using the (fc, r) to denote that space M(ro(l), v^j,, k, A) 
and using w to denote bijectively corresponding spaces we have: 

• Since all v^^l are 24th-roots of unity we have trivially: (fc, r + 12) = {k, r). 

• A composition of Maass operators which raises or lowers the weight by 12 will preserve 
the multiplier system. Hence (fc + 12, r) w {k, r) and we may assume that k € [0, 12]. 

• K is a bijection from (fc, r) to (— fc, — r) — (— fc, 12 — r) w (12 — fc, 12 — r) so there is no 
restriction to assume k G [0, 6] and r g {0, 2, 4, 6, 8, 10}. 

• By using the raising operator we see that (fc,0) w (fc + 2,— 2) w (fc + 2, 10) and by 
repeated use we see that with out loss of generality we can also assume r — 0. 

We are thus justified in our choice of restricting the investigation to the spaces M(ro(l), k, A) 
for k € [0, 6] and v = v^'^l = v'^'^ . 

4.5. Hecke operators for non-trivial multiplier systems 

We know that Hecke operators play an important role in the understanding of the theory of 
both modular forms and Maass waveforms at integer weights. 

For general real weights, the Hecke operators are not important (and may well be non-existent), 
but we will begin with a general definition anyway, and then we will consider two special cases: 
ro(4) with the theta multiplier system and weight ^ and ro(l) with the eta multiplier system 
and weight 1. 

The general introductory discussion is based on [67] but the specific example of integer weights 
(as worked out in detail in [66]) is based on ideas from [76] and [70j . 

Let r C P5i(2,R) be cofinite and suppose v : is a multiplier system of weight fc e M. 

The commensurator of F, comm (F) , in PSL{2, R) is defined as 

comm (F) = jo; e PSL{2, R) \ aTa^^ n F has finite index in F and aFa"^} . 

We know that the Hecke operators are associated with the members of the commensurator, or 
actually with the double cosets, FaF, for a G comm(F). Fix a G comm (F) . It can be shown 
that if the multiplier system v satisfies 

(4.8) v{g) =v°'{g),yg eTna-'^Ta, 

then we can define an associated multiplier system, Va, on the double coset: 

Va : FaF , 

by setting 

Va (91^92) = CTfe {gia, g^) Ok {gi,a) v{gi)v (^2) , 

for all 51, ,92 G F. It might be the case that there exists an associated multiplier system of = x^, 
where x is a character on FaF, even though there does not exist an associated multiplier system 
of V itself. Suppose that Va exists as above, and that we have FaF = uf^^Fai. We then define 
the Hecke operator : M (F, v, fc. A) — > M (F, v, fc. A) by 

d 
i=l 

When F = Tq{N) (or any congruence subgroup of level N), one usually constructs Hecke oper- 
ators T^^ corresponding to positive integers n in which case a — {q^) (cf. e.g. [3], [HI ch. 5], 
[55} ch. 9], [44] or [59] for more details). 
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4.5.1. Hecke Operators for the Theta multiplier System. Consider the case T = ro(4), k = ^ and 
V ~ V0. Let n be a positive integer and a ~ (q n) ■ Then g S ro(4) n a~^To (4)a can be written 
5 ~ ( c/n 'd") I ^iid aga~^ = ( " ^) , with ad — be = 1 and c = mod 4n. It is easy to verify that 

v°'{g) = v{aga~^) and v (g) = v{aga~^) (^) ^ , and hence v'^{g) = v{g) if and only if (^) = 1. 
By (|4.8p this relation must hold for all odd integers d, and hence it is clear that the extension Va 
exists if and only if n is a perfect square. It is also easy to verify that in this case the multipHer 
system is given by {giag2) = u (gi) w (32)- 

Suppose now for simplicity that n = p^, with p ^ 2 a prime number. The + p different 
coset representatives of F in TaT are given by = ^ q ^2 ^ , 6 = 0, . . . — 1, fi = (^01 

and (3h — (op) ' ^ ^ l,...,p— 1. By factorization of these representatives we see that Va (c) 

Va (ctb) = 1, & = 0, . . . — 1 and Va iPb) = for b ~ 1, ... ,p — 1. Hence, for p ^ 2 

( v^-l p-l 

(4-9) - ^(E/(^)+/(^^^)+^E 

For p = 2 the 4 coset representatives are given by = (04),^ = 0, . . . , 3 and we obtain 
precisely the operator in l|4.2p . The above construction is analogous to [6OI pp. 450-451, thm. 
1.7]. Suppose that / € M (ro(4),we, \, A) has the following Fourier series at infinity 

(4.10) !{z) ^ 4=W^i.9«(n),«(47r|n|y)e(nx) 

and that T"" 1/(2;) has a similar Fourier expansion but with coefficients b^P'>{n). Using the for- 

p 1 2 

mula for the standard Gauss sum, J2b=o^ ^ (p^) ~ if else 0, and the twisted version 

EC! (^) e (f ) = (?) (cf- [3 Satz 7, p. 375] or [39l pp. 83-87]) we get that for all 
non-zero integers n 




b^P^n) = |a(p2n) + a(^^^ Qa(n)|,p^2, 
6(2) (n) = a (in). 



and 



(We use the standard convention that a (a;) = if a; ^ Z.) It is now obvious that our Hecke 

operator T"" 1 is equal to the corresponding Hecke operator defined in [291 p. 199]. Observe that 

p ' 2 

our Fourier coefficients a{n) =v^xKatok-Sarnak's Fourier coefficients b{n). 

As usual, we consider Hecke eigenforms in M (To{4:),vg, ^, A) which are simultaneous eigen- 

functions of all T'"^ 1 with p ^ 2. An additional commuting normal operator can be chosen as 

p ' 2 

either or L (these two operators does not commute in general). 

2 

As it turns out, the operator L is particularly useful in connection with the Shimura corre- 
spondence on the Kohnen space, , where L has eigenvalue 1 (cf. Section [6]). 

Observe that the Hecke eigenvalues in this case does not equal to the Fourier coefficients. In 

fact, suppose that / as in (|4.10p is an eigenfunction of all Hecke operators with T^2 if ~ \f 

p ' 2 
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and that a{t) ^ for some integer t. It is then easy to see that 
(4.11) A, = <;^ + M^p^2,and 



A, 



a(t) ^ 
a(4t) 



Using multipHcative relations of the Hecke operators one can prove (cf. |60[ p. 453]) that if t is 
square free, then 

a{tm'^)a(tn^) = a(t)a{tni^n^), for (m, n) = 1. 
Furthermore, if / is also an eigenfunction of then 

2 

a{m)a (4?i) = a(4m)a(n), m, n G Z. 

4.5.2. Hecke operators at integer weights and Fourier coefficients. Consider the modular group, 
r = PSL{2, Z), together with the integer weight fc ^ mod 12 and multiplier v — v'^ . 

It can be shown that for all positive integers n and m with kn = —km = k mod 12 we can 
construct Hecke operators T^^j. and Hecke-Hke operators acting on M(r, w,fc. A). Using 

these operators one can obtain multiplicative relations for Fourier coefficients similar to the 
weight zero case. This is shown in detail in [66] and here we will only state the theorem and give 
a brief outline of the ideas of the proof. 

Theorem 4.7. Let k he an integer k ^ mod 12 and R > then there exist a basis of 
M(r, w, A:, i?) consisting of Maass wave forms f with Fourier expansions at infinity 



k 
12 



k 
12 



where the coefficients c{n) satisfies the following multiplicative relations i/c(0) ^ 0. For positive 

T^. "1 = T arid D = 



integers m,n with 12m — 12n — mod k set mi = ni = and D = tttty. // 



(mi + l,£l) = {ni + l,D) = 1 

(4.12) c{m)c{n) = c(0) ^ Xk (d) c 

0<d|(mi + l,ni + l) 

and if (mi — 1, D) ^ (ni — 1, D) = 1 then 



k /(mi + 1) (ni + 1) 



12 V rf2 



(4.13) ci-m)c{-n) ^ AkMciO) J2 Xk{d) c 

0<(i|(mi-l,ni-l) 

where Xk {d) = i'^i'^^^^ (— (^)'^ for odd d) and 



k / (mi — 1) (ni — 1) 



12 V d2 



(4.14) Ak,B^ 



n-=i(j(j-l) + i + fi')'. k = 2l, 



In particular, we see that if A:|12 we have D — 1 and the multiplicative relations l|4.12p and 
l|4.13p are vaHd for all positive integers. 

Remark 4.8. As in the weight zero case and the coefficient c(l) one can show that if £) = 1 
and / is an eigenfunction of all Hecke operators (defined below) then c (0) ^ unless / (z) is 
identically 0. In the case Z? > 0, if c (0) = and / (z) ^ we can choose an integer no such that 
c ('^o) 7^ and obtain a similar set of multiplicative relations. 
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The proof of the "positive part" of the theorem, i.e. equation (|4.12|) reUes on the construction 
of a family of Hecke operators j. with km = k mod 12. It is shown that this family consist of 
self-adjoint operators commuting with each other and the weight k Laplacian. An explicit form 



of T^.k is 



az + b 



(4.15) T^J(^)= 1 J2 Xk{d)J2v{T)''f 

^ ad=m,d>0 6=0 

It is not hard to show that if / {z) has Fourier coefficients c (n) then T^^^ has coefficients 

, / nm kim — (P) \ 



h{n)= J2 

0<d\ ^m,n- 



12^2 



from which we see that if j./ = A,„/ and c (0) 7^ then 



An 



c(0) 



k (m - 1) 
12 



Xk {D)c 



k (m - 
12D2 



The proof of l|4.12p is concluded with a proof of the following multiplication rule (using straight- 
forward calculations and induction). If km = kn= k mod 12 then 



rpV rpV 



E Xkid)T^^ 

0<d\{m,n) 



To obtain the "negative part" of theorem, i.e. I|4.13p we have to consider another family 
of operators, 6J'^ j. with km = —k mod 12. These operators are given as a combination of a 
bijection Q ~ J mapping M(r, w, fc. A) to M(r, zJ, k, A) and an Hecke operator T^^ ^ mapping 
M(r, w, k, A) back to M(r, v, k, A). Here E'^ = E^^^ o ■ ■■ o Ef^^^ o E,: maps M(r, v, k, A) to 
M(r, w, — fe. A) bijectively since \ > \ and J reflects this space back to M(r, v, k, A). is given 
by l|4.15p with v interchanged with v. The operator Q is similar to the operator defined by Maass 
in p. 181]. Using ^^^-^J^ it is easy to show that O^f = Ak,Rf for all / G M{r,k,v,X) 
and that if / (z) has Fourier coefficients c{n) then j,f (z) has coefficients 



d{n) 



E 



Xk id) d 



, I nm k [m + rf^) 



d2 



12d2 



where c' (n) — c {—n) 



1, 



0<d\ ym,n— 

n > 1, 



Afe,_R, n<0. 



If J = Mm/ and c (0) ^ then 



:(0) 



k [m + 1) 
12 



Xk {D)c 



' k{m + D'^ 
12D2 



The multiplication rule for the operators 0^ is that for km = kn = k mod 12 we have 

0<(i|(m,n) 



This in combination with the expressions for fim and Am concludes the proof of l|4.13p . 
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Example 4.9. Look at the specific case k = 1 and a function / G M(r, u,l,A) then by by 
setting m = 1 in l|4.12p and using the normalization c(l) = 1 we see that 



c{n)^ X! Xi{d)c 

0<ti|(12ri+l,13) 



(I3(12n+ 1) - (f 



l2dP I 

and hence if (12n + 1, 13) = 1 then we get a striking proportionahty relation: 



(4.17) c(n)=c(0) c(13?i + l) + c 



(4.16) c (n) =c(0)c J^Mi^!L±l) 1 j ^ c(0)c (13n + 1) . 

For 12?i + 1 = 13^ we get 

Is 

and if (/, 13) = 1 we can combine these two equations and obtain: 

(4.18) c (n) =c(0) (c (13n + 1) + c(0)c (n)) 
and hence 

(4.19) c(n) = ^^^c(13n+l). 

Consider now m = -1 in (|4.13p and note that Ai^r = -i?^. If (127i - 1, 11) = 1 then 

(4.20) c(-n) = -^^c(lln-l), 

c(~l) 



and if (12n - 1, 11) = 11 then 
(4.21) c(-n) 



c(-l) 



c(lln- 1) 



11 



We have not seen relations of type l|4.16p - l|4.2ip earlier in the literature. For some numerical 
examples see Tables [2] and [3l Examples of relations for higher weights can be found in [66] . 

Remark 4.10. An alternative approach to the above coefficient relations in the case of weight 1 
is to identify M(ro(l), w^, 1, A) with a subspace of M(ro(144), (^) ,1,A) via the map /(z) i-^ 
g{z) = /(12z). Cf. e.g. [511 sec. 2.4.7]. 

This identification also provides an explanation for the occurrence of CM-type forms found 
(numerically) in M(ro(l), u^, 1, A). These forms have eigenvalues in an arithmetic progression: 
Rk = ^T^^Tify, for fc e Z. See Table H 

5. The Eisenstein series for PSL{2, Z) at weight zero 

In case there is a cusp pj at which the multiplier system is singular (i.e. v{Sj) ~ 1) we have 
a continuous spectrum: [j, oo) (with multiplicity equal to the number of singular cusps), and in 
general we can not say much about the embedded discrete spectrum in [|;,oo). 

Examples of singular cusps are the cusp at infinity for the eta multiplier and weight k = 
mod 12 on PSL{2, Z) and the cusps at and ioo for the theta multiplier and weight ^ on 

ro (4). 

Remember that Maass waveforms are part of the discrete spectrum, but as we continuously 
"turn off" the multiplier system, i.e. for v = w^*^ we let fc ^ 0, the continuous spectrum will 
emerge in the limit. For this reason we want to review some details concerning the Eisenstein 
series on the modular group. 
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At weight and singular character X; the continuous spectrum of A = Aq is the interval 
[j, oo) and the eigenpacket is given by the Eisenstein series E{z; s; x) defined by 

E{z;s;x)= E x{T-')mTz)r, 
Ter^\r 

where Too ~ [S]. For the trivial character we have the Fourier series expansion (cf. p2 p. 65 and 
p. 76]) 

E{z; s; x) = + f{s)y^^' ipn{s)^/yK^_i{2TT\n\y)e{nx), 



where 



fis) = ' ' . , and 

r(s) C(2s) 



r{s) c(2.s) ■ 

Hence we can see that for s = ^ + iR , the nth Fourier coefficient of E{z; s) is given by 
(5.1) c(n) = ifn(l + iR) =K ■\n\''^a^2^R{\n\) 



2 

d||n|,d>0 

where K = K{R) is a constant dependent on R. For a prime p we get 
c{p) = K ■ p'"{l + p-"^'") ^ 2K ■ cos(i?lnp). 

Using this formula we can compute quotients of various c{p) (e.g. ^||y) and compare this with 
corresponding quotients for the experimentally obtained forms of weight fc « 0. 



6. The Shimura correspondence 

We know that the 0— function is an automorphic form (not a cusp form) of weight ^ on ro(4), 
hence we consider ro(4) together with the 6'— multiplier system (cf. section \2^ . 
6.1. Introduction the holomorphic case 

We will consider the Shimura correspondence only in the particular case of trivial character and 
level a square-free multiple of 4. Let Sk{N) denote the space of holomorphic cusp forms of weight 
fc e Z (and trivial multiplier) on ro(iV), and let Sf._^_x{AN) denote the space of holomorphic cusp 
forms of weight fc + i, fc e Z and multiplier vg, on ro(4iV). The Shimura correspondence is a 
correspondence between the space S'^+i (4-/V) and spaces S2k{N') for certain integers A^'|4A^ (e.g. 
TV' = 2N or N). 

The map from Sf.^i to S2k was first constructed by Shimura [60] and later an adjoint map 
from S2k to Sk+L was constructed by Shintani (62] ■ The former uses a Dirichlet-series and 
the latter uses an integral against a theta-function. Both these maps commute with the Hecke 
operators that are acting on S2k{N) and S'^.+i (4A^) respectively. Kohnen [32l[33] proved that for 
N odd and square-free, the correspondence is a bijection between the newforms on S2kiN) and a 
certain subspace, C Sk+i (4iV). The subspace is composed of Hecke eigenfunctions whose 
Fourier coefficients, c(n), satisfy certain vanishing properties; namely, c(n) = for n = 2, 3 mod 4 
(see also Section fO]) . 
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Following from the Shimura correspondence is also a connection between certain Fourier coeffi- 
cients of the half integral weight forms and critical values of twisted L-series for the corresponding 
integral weight form. Cf. e.g. [73 1 173]. [SS l ISi] and [61]. 
6.2. The Shimura correspondence for Maass waveforms 

The extension of the Shimura correspondence and Kohnen's result to spaces of Maass wave- 
forms has been investigated by e.g. Sarnak [57], Hejhal [2^, Duke [14], Katok-Sarnak [29j . 
Khuri-Makdisi [S^, Kojima [33 [371 [M] , Biro [^ and Arakawa [2]. Of these, [30] and [37] [38] 
are written in the more general setting of Hilbert modular forms for number fields. Reading [30] 
together with [61] and [29] gives a good picture of the current state of affairs. 

Throughout this section let M {N, R)~^ denote the space of even (with respect to J : z i-^ 
—z) Hecke normalized Maass waveforms in M (ro(A^), 1, 0, i?) and let Mi(4, i?) denote the 
space of Hecke normalized (with respect to all Tp2, p ^ 2) weight ^ Maass waveforms in 
M (ro(4), we, i, i?). Also let V+ C (4, R) be defined as in Section [431 

For ./ e M {N, and (/) e Mi (4, R) we will use the following notation: 

OO 

/(z)= ^ A{n)^K,R{2TT\n\y)e{nx), 

n— — oo 

and 

0(2)= Y] -jf^Wi (^^^ {47r\n\y) e (nx) . 
\/\n\ ■* " ^ ^' 

n— — 00 V I I 

(Both expansions are given with respect to the cusp at 00). 

The existence of a Shimura correspondence and an inverse for Maass waveforms is expressed 
by the following proposition (essentially [30l thm. 5.1 and 5.2]): 

Proposition 6.1. 

a) The Shimura correspondence gives a map $ : Mi (4, i?) ^ M (2, 2R)^ . 

h) Conversely, let f G M (2, 2i?)+. Then there exists a (f> e e Mi (4, R). 

c) If A(p) is the Hecke eigenvalue of f with respect to the operator Tp and Xp is the 

eigenvalue of 4> € $~^(/) with respect to T'"^ 1 then we actually have 

p > 2 

A{p) = Xp. 

Proof. See the proof of [301 thm. 5.1 and 5.2]. For c) observe the difference in normalization of 
the Hecke operators. Cf. also [36l thm. 2], Proposition 16.31 a) and l|4.1ip . □ 

Remark 6.2. From [29l prop. 4.1 and 2.3] we know that the above correspondence $ actually 
restricts to a map between M (1, 2R)'^ and the subspace . And hence a new form in M (2, 2R)~^ 
will be mapped to V~ , the orthogonal complement of V'^ . Vice versa, a form in V~ will be 
mapped to a new form in M (2, 2i?)^. 

To generalize the results mentioned at the end of the previous subsection to Maass waveforms 
we need the following definition. For / e M(2, i?)+ with Fourier coefficients {A{n)} and a given 
Dirichlet character x we define the x-twisted L-series of / by 

00 

i(/,X,s)-^^(n)x(n)|n|-^-i 
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Proposition 6.3. Let (f> G Mi (4, R) have Fourier coefficients {a{n)} and correspond (via Prop. 

Tip to f = $(0) e M(A^,2i?)+ (where N = 1 if and only if (f) £ V+ , otherwise N = 2) with 
Fourier coefficients {A{n)}. Let t G Z+ be square free and let x't be the quadratic residue symbol 
(-) considered mod Nt (i.e. we have X(("-) = («) */ ("j ^) = 1 '^^'^ Xt("-) = otherwise. 
Then the following properties hold: 

a) If o-ii) 7^ then A{n) can be expressed by: 



1 2 \ 



(6.1) A{n) = 2^ '^'y n e 



kd—n 

fc>0 



x't (fc) « {i<P 

Vk a{t) 



To express a{t) in terms of Ain) 's we get three cases. 

b) If f is an oldform. i.e. / e M(1,2_R)+ then 4> G and hence 

a{t) = 0, i = 2,3 mod 4. 

c) If f is a newform with eigenvalue e with respect to the involution ^ then 

a{t) = 

for t = 5 mod 8 if e = I and for t = I mod 8 i/ e = — 1 . 

d) For all other (square free) values oft the following formula holds 

\aitf^Q^Lif,xuO) 

where Q is a constant independent oft. 

e) If 4> is a normalized Hecke eigenform (for all T'"^ i , p an odd prime) then f is also 

P , 2 

a normalized Hecke eigenform. 

Proof Cf. [601 PP- 448(prop. 1.3), 458 (main theorem), 474 (line -11)], and [30l p. 422] for the 
choice of x't (the N = 1 case is simply to incorporate the results of [29]). Relation a) follows 
from |30l thm. 5.1 (2)], which in our case can be written as 

(6.2) J2 ^(")""' = c • L ( . + -, xi) 5] a n-^ 



=1 



(Cf. e.g. also [601 P- 458], ^ p. 159], [SH prop. 3.1].) The relation b) follows from |29l prop. 
2.3] (see also |32]) and relation c) is implicit in [30l thm. 8.1 (8.6)]. Simply observe that in our 
normalization the relevant factor in the product is given by 



(V2A{^. 



v--'^'-e))-(-(l 



since ^(2) = for a newform (analogous to [3l thm. 3]). (Alternatively consider the sign of the 
functional equation for L (/, xt, s).) Relation d) is the Maass waveform-analogue of [35l thm. 1] 
and follows from |30[ thm. 8.1] (in the notation of |30j we have a{n) = \n\i^{n; </>)). □ 

Remark 6.4. That the map $ is well-defined and surjective from Mi (4, R) onto M (2, 2R)^ 

follows from (301 thm. 5.1 and 5.2] and "multiplicity one" for the Hecke operators Tp on M (2, 2R)^ 
(cf. [IHl thm. 4.6]). Experimentally we have observed that <f> restricted to also seems to be 
injective. Theoretically, this is still an open problem that might be possible to resolve using the 
trace formula for the Hecke operators Tp2 on . 



COMPUTATION OF MAASS WAVEFORMS WITH NON-TRIVIAL MULTIPLIER SYSTEMS 



19 



Remark 6.5. Note that the same argument as for 16.31 c) implies that to any / G M(2,2_R)''' 
with Fourier coefficient A{2) = ^ there corresponds a function (j) € (/) C Mi (4, R) with 
coefficients a{n) = for either n = 1 or 5 mod 8 respectively. Hence, since oldforms on ro(2) 
occur in pairs, we can choose two forms /i,/2 with Ai{2) = ^ respectively ^42(2) = These 
two functions thus correspond to two Hnearly independent non-zero functions in Mi (4, R) which 
hence is at least two dimensional when 2R is an even eigenvalue for PSL{2, Z). 

The relations a)-c) in Proposition 16.31 appeared in |64[ ch. 2] as experimental observations 
(cf.[83 in particular Tables [MH]). 

See also [H §4.1] and O p. 633], and (HH pp. 502 (bottom) - 504 (top)] for some additional 
perspectives. 

7. Some computational remarks 

We recall the key ingredients in the standard Hejhal's algorithm (cf. e.g. [HI [261 [U] , [65| [63] 
and [6]) to compute weight zero Maass waveforms on cofinite Fuchsian groups. First of all the 
asymptotic properties of the K-Bessel function are used to obtain an approximation the Maass 
waveform by a truncated Fourier series. By viewing this as a discrete Fourier transform one 
can use the inverse transform to express any coefficient as a linear combination of the other 
coefficients. Finally, by using the assumed automorphy properties of the function one obtains a 
non-trivial linear system that is satisfied by the Fourier coefficients. 

To use this algorithm to also locate eigenvalues the most general method is to use two different 
sets of sampling points for the inverse transform and try to minimize the difference between the 
correspondingly computed solution vectors. 

The following two modifications are needed in order to make the weight zero algorithm work 
in the general case: 

• the i^-Bessel function needs to be replaced with the Whittaker function, and 

• the automorphy condition needs to incorporate the multiplier system. 

The first modification, although trivial in theory is the most complex in terms of the numerics 
involved. There was no efficient algorithm for the Whittaker function in the literature and a new 
algorithm had to be developed. We used the integral representation (cf. [22l p. 375 (top)]) 








— fc; -; .t(1 + cosht) 



cosh [iRt] dt, 



where ^' is a confiuent hypergeometric function together with a stationary phase method to 
develop a robust and efficient algorithm. This algorithm is in essence similar to Hejhal's algorithm 
(cf. [25]) for the K-Bessel function, Kift{x), and the generaHzation of it made by AveHn (cf. [U 
AveHn]) to a complex argument, Ks[x), s £ C. The details of this algorithm is described in [SU 
Ch. 4]. 

Let F be a Fuchsian group with fundamental domain J- and let pj, aj,\ < j < k, qi and 
Ui, 1 < i < kq he as in Section [L2l Let I (w) = i if qi is the closest (in the hyperbolic metric) 
parabolic vertex to w G J-. 

Consider z e 3i and let Tj G F be the pull-back map of ajz, i.e. TjUjZ G T. Observe that 
/j(^) = /|<Tj {z) = j'y, k)~^f{(Jjz), and with z* = aJ^^-^Ui(j)TjajZ where we set / (j) = / (Tjajz) 
(cf. |65[ p. 23]) the automorphy condition now becomes: 

f,{z) = j,^{z;k)-^f{a,z)^j,^{z;k)-^f{T-^U-]^aj^,)Z*) 
(7.1) = ja^{z; k)-^v{T^^UYf^^yk)jj,-ijj-i^ k)f{ai(^j)Z*) 



Hi) 



20 



FREDRIK STROMBERG 



The entire setup of the Phase 1 algorithm, i.e. locating eigenvalues and computing the first few 
Fourier coefficients, goes through exactly as in the case of weight (cf. [65] or [64} Ch. 1]) with 
some trivial modifications. For the sake of completeness we will give an outHne of the modified 
algorithm. Consider / G M(r, v, k, A) and using the notation Ui = n + Ui f has a Fourier series 
expansion at the cusp pi : 

°° c (n) 



— OO VI ^1 

and since the Whittaker function is ultimately exponentially decaying, given an e > 0, there 
exists a constant (depending on y and e), M{y), such that 

/,(z)=/,(z) + [[£]], 

where we use [[e]] to denote a constant of magnitude less than e. The truncated Fourier series 

Ci[n) 



~M{Y) 

is now viewed as a discrete Fourier transform, and if we take the inverse transform over the 
horocyclic points: z,n = Xm + iY , 1 — Q < m < Q, where Xm = 2^(5 ~ ttt-) we get: 

I Q ^ I Q ^^^^^ c-n) 

'^^ l^Q m=l-Q ;=-M(y) V K^l 

i=_M(Y) Vl'^l 1-Q 
C' ( tl] 

=^W^|.,„(„.)..fi(47r|n.|r), 



where we used the fact that 

Now we also have fi{zjn) = Xmili{jn,i){z*mi)^ where (by {HH)) 

Xrm = jaAZm] C/-^,, <t,(,„,,, (^™; k)w{T~^UY^^^^yai^,^^,))v(T~^UY^^,^yk). 

Hence 

c (n) 1 



2Q 



1 ^ 
77 Yj 



(m,i)\Zjni)G[ TiiXm) 



K Mo 

= Y Y c,(0K^+2[[e]], 



j=l l=-Mo 
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where 



1 1 



Q 




We then define 



;f by 



1 



and if we neglect the error of magnitude e we finally obtain a (well-conditioned) linear system 



which must be satisfied by the Fourier coefficients of /. Here V is the matrix V^j and C is the 
vector Ci{n), both depending on R and Y. The basic idea of Phase 1 is now to locate eigenvalues 
R together with sets of corresponding Fourier coefficients, Ci{n), by solving (*) repeatedly for 
different R's, and seeking those values of R for which the solution vectors C = C{R, Y) are 
stable under changes of Y. That is, if C{R,Yi) sa C(i?, I2) (for Yiand Y2 < Yq, for some 
suitable constant Yq) we take it as an indication that the R is close to an eigenvalue and that 
the components of C are close to the corresponding Fourier coefficients. For more details and 
justifications see [Ml sect. 3.2]. 

We note that the Phase 2 algorithm (i.e. computation of a larger set of Fourier coefficients) 
also generalizes to non-zero weight in a similar manner (cf. |65[ sect. 3.3]). 



The experimental excursions have been directed towards three essentially different subjects, 
but, in each, we have worked in an exploratory spirit. 

• First we tried to get an over-all picture of the distribution of small to middle-range sized 
eigenvalues on P5i(2,Z) (and the eta multiplier) for "large" weights, e.g. R G [0,14] 
and k € [0.1,6]. 

• Second, we continuously "turned off" the multiplier system vf^^ on P5i(2,Z) by letting 
the weight A: ^ and studied the varying distribution of eigenvalues, Nk{T), as well as 
the formation of the continuous spectrum. 

• There are some cases where arithmeticity plays a role even in nonzero weight. We studied 
the Shimura correspondence between weight zero forms on ro(2) and weight one half 
forms on ro(4). And we also studied weight one forms on PSL{2,'Z) which correspond 
to Hecke eigenforms on ro(144) with a Dirichlet character. 

8.1. Varying weight 

The first experiment considered was to tabulate the first few eigenvalues (up to i? = 14) for 
PSL{2,Z) and the multiplier system u^'', of weight k £ (0,6] (cf. Section BXTj) . We made the 
computations for k £ [0.1, 6] with a grid size of 0.01, and the results are presented in Figure [TJ 
We stress here, that the arcs in Figure [T] terminate at /c = 0.1; it is not excluded (and actually 
expected) that they might go lower. 

For some examples of eigenvalues for "large" weights see Table[TJ This data should be compared 
with data obtained by Miihlenbruch, [46l p. 160], who used a completely different method (with 
much less accuracy). We note here that as R increases, the negative Fourier coefficients seem to 
grow rapidly in magnitude (as compared to the positive ones, with the normaHzation c(l) = 1) 
for "large" weights. 

We believe that we have found all eigenvalues with (i?, k) G [0, 14] x [0.1, 6]. This belief is 
founded on the "continuity" of the resulting graphs Rj{k) (cf. Figure [T|), where Rj{k) is the j-th 



(*) 



8. Numerical results 
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Figure 1: Section of eigenvalues with < i? < 14, and weight 0.1 < fc < 6. 



OS 

6 




1 2 3 4 5 6 

Weight *: 




eigenvalue at weight k, considered as a function of k. By standard results (e.g. [TOi p. 149]) Rj{k) 
should be a real analytic function in this interval. 

Remember that, for fc > 0, the smallest eigenvalue, Ami„, corresponds to the function F{z) = 
y^ri{z)^'^. A lower bound for the second smallest eigenvalue, Ao(fc), is discussed in [101 P- 183]. 
Bruggeman finds two such bounds, both positive, which he calls /lo(fc) and ^i(fc) {ni is better 
than ^0 in a certain interval / C [0,2].) Figure [2] shows a comparison between the i?- values 
corresponding to these bounds and the smallest experimentally found eigenvalues in the interval 
fc S [0.1, 6]; we see that Bruggeman's bounds can be improved quite a bit. 
8.2. Small weights 

The investigation of eigenvalues for small weights k has been done in the interval R G [0, 20], 
and k e {lO^-' [1 < j < 12} C [lO^^^,0.l] . We beheve that most cusp forms were found. Let us 
use the notation 

A,(fc) 
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for the j-th discrete eigenvalue at weight fc, and 4>j{k) for the corresponding cusp form. It is then 
a basic fact that Xj{k) depends continuously on k, but it can also be shown that for k S (0, 12) 
Xj is even real analytic in k. That is, 0j(fc) belong to an "analytic family" in the terminology of 
Bruggeman (see [lOlIII])- In connection with this, it should also be noted that our experiments 
support the statement in Observation 71 in [46] that the first few cusp forms at weight do not 
belong to an analytic family of cusp forms defined in the interval (—12, 12). Indeed we find that 
in the range considered we actually seem to have Xj{k) — *■ as fc — > which is consistent with 
Bruggeman's result (cf. case ii) b) in Prop. 2.17 in [TOl p. 149]). 

Our experiments indicate that for fixed small k, the "generic" cusp forms 4'j{k) can be divided 
into two classes: 

• C{k), and 

• E{k). 

Here C{k) consists of functions (j)j{k) such that Aj(fc) is close to an eigenvalue of a cusp form at 
weight and the Fourier coefficients are close to the corresponding coefficients of the weight-zero 
cusp form. 

E{k) on the other hand consists of functions (j>j (fc) such that the Fourier coefficients are close 
to the Fourier coefficients of the Eisenstein series E{z, s) where Xj{k) = s(l — s). 

The typical difference between the Fourier coefficients at weight k and weight are in both 
cases basically of order k; for the forms in C(fc), the distance between Aj(0) and the corresponding 
discrete eigenvalue at weight is much smaller than k. 

The "generic" in the statement above means that we exclude certain places where the families 
0j(fc) changes character between C(fc) and E(k). In these cases we have a situation of almost 
multiplicity 2, and in too coarse resolution it actually looks Hke the analytic famihes intersect. 

Weyl's law For non-trivial 77-multiplier and a fixed non-zero weight k £ (0, 12) on PSL{2, Z) 

is 

(8.1) Nk{T) = < T, weight k} = ^ - -In 1 - e"^' + SiT) + 0(1), 

for T > 1 (cf. [2II p. 466] with kq = 0). Our experiments seem to suggest that as fc ^ 0, the 



and indeed it is easy show that the 



main contribution is proportional to the factor In 

0(1) term is even uniformly bounded in k. 

To obtain asymptotics for S{T) when A: ^ (and T is bounded) is a bit more involved. We 
want to generalize \22\ thm. 2.29, p. 468] by following the pointers provided in |221 proof of thm. 
2.29, p. 468]. Basically, our goal is to single out any terms that grow as A; in the estimates 
for S{T). 



By careful bookkeeping we see that these terms are all of the form In 



1 - e" 



and this kind 



of terms only show up when we apply the functional equation for the logarithmic derivative of the 
Selberg zeta function, Z (s) ( |22[ thm. 2.18, p. 441]). To obtain the necessary estimate of ^ (s) 
we study the integral in the left hand side of \2T\ prop. 8.6, p. 121]. Let the assumptions of \21\ 
ass. 8.5 p. 120] and |2H p. 125 row 5] apply, e.g. s = a + iT, U = r + 10 and ^<^ + e^(Ji<a. 

Through the integral estimates \21\ prop. 8.7-8, 8.10-11] we consider either ^(^) — 2 In 1 — 

or ^ (^) , depending on whether the functional equation is appHed along this part of the path 
or not and then collect the extra terms at the end. It turns out that we consider the modified 
integrand in all parts to the left of 5R(^) = ^ — e. With this modification, all estimates in 
\21\ prop. 8.7-8, 8.10-14, thm. 8.15 and 8.17] go through except for the addition of a term 
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Rniis) 



In X 



e e 



/m(s) where /„ 

Imis) = 



s) is the integral 



1 
27ri 



7 - 

where the path is 7 = 7^; = - e - ill, - A- iU]Ll[-A - iU, - A + iU]u[-A + iU,^ - e + ill] . 
The integrand has no poles to the right of 3? (^) = 5 — e and is hence equal to the integral from 



ill to 



\Imis)\ < 



+ ill. Elementary estimates now show that 



27r 



-u £2 + (t - ry 



rdt < 



-2e 



-u 1 + 



rdt < 



Using X 
see that 



Ti (cf. [211 p. 138, row -6]) together with e = (cf. [HI ass. 8.16(b), p. 130]) we 



\RrrL (s)| < a:: In 



1 



In; 



In a; 



■In 



1 — e 8 



This error term is now to be added to in |2H thm. 8.15] and subtracted from the right hand 
side of [m thm. 8.17] and finally in [HI p. 135 (8.18)] we need to add - /? S [R„i (s)] to the 
expression of nS{T). Thus the analogue of [HI thm. 8.1, p. 119] and [HI thm. 2.29, p. 468] 
is S (T) ~ O (1^) + 0(1) In 1 — e^* (with implied constants independent of k) and inserted 
into l|8.ip we conclude the following Weyl's law: 



i.2) 



Nk(T) 



12 



In 



e 6 



T 



Oil) 



O 



T 



uniformly for < k < 0.1 (say) and T > 2 and the 0(1) term is bounded in magnitude by 
^e"^ w 0.065. 

We have computed Nk{T) experimentahy for k = lO^-'jj = 4, . . . , 12, and < T < 20. Figure 
[3] shows a picture of the experimental values compared to the theoretical values obtained by l|8.2p 
(with the O-terms neglected) and the difference indeed seems to be constant over this range. 

From the form of the Weyl's law above, we also see that the successive spacings A„(A:) = 
Rn+i{k) — Rn{k) should look about like -dw^, i-e. A„(fc) w — — yj\ — kTi ~ 1 , \^ 1 as /c 0. Figure 

[4] provides a nice illustration of this fact, where it is clearly seen that the average spacings are 
almost constant for small k and this constant is roughly proportional to ji^- 

The mean gap over 100 typical cases (i?„, 5 < n < 155) turned out to be 0.117 at /c = 10^^, 
0.107 at k = 10^1", 0.101 at A: = 10"" and 0.095 at fc 10"^^ The relative quotients, 
^7i{k) [^ + ^ \^^^\] ' are 1.019, 0.990, 1.001 and 0.9934, respectively and it is not inconceivable 
that one obtains 1 in the (logarithmic) limit. 

8.2.1. Level repulsion. From figures like|l]one may be tempted to think that there are horizontal 
hues corresponding to cusp forms at weight which crosses the hues that are going down (i.e. 
corresponding to the Eisenstein series at weight 0). This is not the case! If we look closer we will 
see that there is actually "level repulsion" here, i.e. the horizontal "cusp-form-line" turns down 
before the "near crossing" and becomes an "Eisenstein-series-line" and the previous "Eisenstein- 
series-hne" turns into a "cusp-form-hne". See also Figures [5] and [6l More precisely formulated: if 
there is a "near crossing" at the weight fco close to the eigenvalue i?o ~ Rji^o) ~ Rj+i{ko), then 
there are two analytic families </>j(fc) and (l)j+i{k) such that for some S > e > 0: 

, l'c(fc), ke[ko + e,ko + S], 

[Eik), k£ [fco - (5, fco -e], 
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Figure 3: Plot of Weyl's law with constant T = 20 and weight fc — > 
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Figure 4: Section of eigenvalues with 9 < i? < 14, and IE ~ 9 < k < IE - 7. 
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and 

, (E{k), ke[ko + e,ko + 5], 

^'^^^ \C{k), fee [fco-<5,fco-e], 

and in the interval (fco — e, /cg + e) both families display a mixing between the two types E{k) 
and C{k). In fact, the Fourier coefficients of (/>j+i converge (as fc ^ fco — e) to values close to the 
Fourier coefficients of 4>j for fc > fco + e and vice versa. Since the two functions also need to be 
orthogonal it is clear that the Fourier coefficients exhibit "wild" behavior in the small interval 
surrounding the "near crossing". Note also that, as fc ^ 0, all <i)j{k) e E{k) and Rj{k) 0. 

See Table [9] for examples of Fourier coefficients corresponding to eigenfunctions of types E{k) 
and C(fc), close to an avoided crossing at weight fc = 9.044605824E' — 8. Table [TOl illustrates 
the agreement between the Fourier coefficients of a more generic cusp form in E(k) and the 
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Figure 5: Level repulsion at R = 13.779 



14.0 - 

13.9 ■ 

13.8 ■ 

13.7 - 

13.6 - 



l.Oc-03 1.2C-03 1.4C-03 1.6c-03 1.8l--03 
Wciahl A- 



1.6360C-03 1.6366C-03 1.6372^-03 1.6378c -03 

Weiahl f; 



1.636794C-03 1.6368OOe-03 1.636806e-03 
Weight f; 



1.60C-03 1.63C-03 1.660-03 

Wcishl A- 



I.63670C-03 1.63676e-03 1.63682C-03 1.63688t:-03 

Weight A- 



1.6367962e-03 1.6367970e-03 1.6367978e-03 

Weight A- 



corresponding coefHcients of the Eisenstein series (appropriately normalized) at weight 0. The 
level of agreement is striking to put it mildly; likewise in Table[9]for the C{k) eigenfunction. The 
"1 for 1" nature of this convergence in the presence of a limiting continuous spectrum seems not 
to have been suspected earlier. Cf. [24^ thm. 6.6 and cor. 6.9] and |23j . 

The fact that the system seems to avoid accidental degeneracies by means of level repulsion 
and avoided crossings is in agreement with the Wigner-von Neumann theorem, cf. [72] . 
8.3. Lifts at weight 1 

As we saw in Section l4.5.2l we could prove the existence of certain Hecke relations at weight 1 
(e.g. (|4.19p . I|4.20p and l4.21]) ). Tables [2] and [3] contain numerical verifications of these relations. 
Table |4] contains a list of computed eigenvalues on M (ro(l), w^, 1, A), and the eigenvalues corre- 
sponding to cosine CM-forms are indicated. In these cases, we have computed the actual error 
since we know the exact eigenvalues: 

27rfc , , 

Rk = — k e Z+. 

In (7 + 2^/12) 

Note that the actual error is in general much smaller than the error-parameter which is 
basically i?(Yi, ^2) = |c(2) - c'(2)| + |c(3) - c'(3)| + |c(4) - c'(4)|, where c(n) is computed with 
Yi and c'{n) with Y2. 
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Figure 6: Level repulsion at R ^ 9.533 . 



9.53366 - 



1.20C-08 1.60e-08 2.00t^-08 2.40i;-08 2.80c-08 
Weighl k 



9.53362 - 
2.4350C-08 



2.4353e-08 2.4356e-0 
Weight k 



t3549i;-08 2.43551c-08 2.43553^-08 2.43555c-08 2.435523e-08 2.435525e-08 2.435526e-08 2.435528c-Q 

Weight k Weight k 

\ \ ~=r 9.533695264 F \ ~ I ^ 



9.53369530 - 



9.533695262 - 



................ ^ ....... . 

1525 - - 



9.53369520 - 



9.533695258 - 



2.4355252e-08 2.4355257^-0! 
Weight k 



2.43552545e-0f 



2.43552547e-08 2.43552549c-0f 
Weight k 



8.4. Half integer weight 

We now consider the case of ro(4) and the ^-multiplier system for weig ht fc = i. The aim of 
our investigation in this case was to study the Shimura lift, and in particular to investigate the 
dimensions of the spaces of half integer weight forms. As remarked at the end of section [6] several 
properties of the Shimura correspondence were observed numerically, and in the original version 
of these notes, [Ml ch. 2], we formulated a number of experimentally inspired conjectures. Except 
for the question of injectivity these conjectures have now been resolved with Propositions 16.11 
and l6.3l (With the obvious correction in the first conjecture: the dimension of !M(ro(4), wg, i, i?) 
corresponding to an oldspace is at least two dimensional.) 

Tables [5] and [6] contain examples of Fourier coefficients at weight ^. Note especially in table [6] 
the agreement with Proposition l6.3l c) displayed by the Fourier coefficients c{n) for n = 5, 1 mod 8 
respectively. Here, it is known that 8.92287648699174 and 12.09299487507860 on ro(2) corre- 
spond to the eigenvalues 1 and —1, respectively, with respect to the involution 0^2 • 

Table [8] contains a comparison of Fourier coefficients computed both from forms on ro(2) via 
l|6.ip and computed directly. Additional Fourier coefficients for the weight forms are available 
in Table [71 
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The normalization we have used here is the usual c(l) = 1. 
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Table 2: Fourier coefRcients for a CM-form 
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^ This quotient is deduced from formula (|4.16p or l|4.17p on 

bc(-l)2 = -i?2c(0)(c(10)-c(0)) 
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Table 3: Fourier coefficients for a non-CM-form 
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This quotient is deduced from formula l|4.16p or 
(|4T7|1 on p.[ll 
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Table 4: Eigenvalues for M{PSL{2, Z), 1, v^) 

R H{yi,y2) True error ^ 

lE-15 ^ 

lE-14 ^ 
4E-13 ^ 

6E-14 ^ 

5E-14 



2, 


,38549209578045 


lE- 


12 


3. 


,66240686698667 


lE- 


■11 


4. 


,77098419156091 


lE- 


13 


5. 


,77698688078694 


6E- 


■12 


6. 


,64285171613711 


lE- 


11 


7. 


,15647628734173 


lE- 


11 


7. 


,82634704540775 


lE- 


13 


8. 


,66620831896793 


7E- 


14 


9. 


,45156176783224 


7E- 


12 


9, 


,54196838312186 


2E- 


13 


10, 


,21802876776059 


3E- 


13 


10, 


,65897262920241 


2E- 


■12 


11, 


,27526358349387 


2E- 


13 


11, 


,92746047890219 


7E- 


■11 


12, 


,15792337422149 


7E- 


13 


12, 


,55403509998720 


lE- 


13 


13, 


,00123950642372 


9E- 


13 


13, 


,67542640643589 


2E- 


13 


13, 


,71353384358095 


lE- 


■12 


14, 


,31295257468268 


lE- 


■12 


14, 


,47039277253940 


lE- 


12 


15, 


,03845367358721 


2E- 


13 


15, 


,39856858348441 


lE- 


12 


15, 


,85705128717333 


3E- 


■10 


16, 


.14536205734475 


3E- 


11 


16, 


,45061260131967 


4E- 


11 


16, 


,69844467046304 


4E- 


12 


16, 


,93043847896222 


4E- 


15 


17, 


,51562192885174 


2E- 


■10 


17, 


,59022138305996 


lE- 


■10 


18, 


,13826107340244 


3E- 


11 


18, 


,32637702205910 


4E- 


11 


18, 


,76341585146817 


lE- 


11 



lE-14 ^ 



lE-13 



^ For CM-forms, the true error is computed 
with respect to the eigenvalue Rk = 

where tjq = 7 + 2^12. 
^ These forms correspond to CM-forms. 
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Table 5: Fourier coefficients for 
/i,2 e M(ro(4),i, 6.889875675945) 

Fourier coefficients for fi G , observe tliat a{n) ~ for n = 2, 3 mod 4 



Ti 




Tl 




A 

4 


U.842197d9d75471 


3 


l.L)18z5z9917145o 


5 


U.loo55ozl4(Jo443 


4 


2. 18968040385979 


8 


0.5d9(J7998524429 


7 


2.06305218095270 


9 


—0.33045049673565 


8 


-1.17157116610978 


12 


—0.41296169213831 


11 


-1.02718719694121 




U.DUiOoOo 1 yooUO < 


1 


007^07^1 ^1 A 

z.zy ( oy < oioi4ooi 


16 


0.30482066289397 


15 


-4.08935474990375 


17 


-0.88689598690620 


16 


3.39248165496032 


20 


0.41418282092059 


19 


-1.31804633824673 


21 


0.50212917023175 


20 


1.87725570455517 


23 


-0.00000000000110 


23 


-2.22265197818114 


24 


-1.04429548341249 


24 


0.58816620330140 


25 


0.28984678984391 








|c(4)c(9) - 


c(36)| = 0.2£; 


-08 



Fourier coefficients for /2 ^ V^^, observe tliat a{n) = for n = l mod 8 



n 


a{n) 


n 


a{n) 





0.000000000000 


19 


-0.936283350934 


1 


0.000000000000 


20 


0.192087703124 


2 


1.000000000000 


21 


1.247834928335 


3 


-0.725665465042 


22 


0.411443268212 


4 


-0.710754741008 


23 


-0.018564690418 


5 


0.456158224759 


24 


-1.297582830232 


6 


-1.835059236817 


25 


0.000000000000 


7 


0.481289183972 


26 


-0.610656711780 


8 


0.707106781187 


27 


-0.179166638197 


9 


0.000000000000 


28 


0.340322845698 


10 


-0.651049038069 


29 


0.001563849679 


11 


-0.470914040036 


30 


1.462745418892 


12 


-0.513122971204 


31 


0.400763051265 


13 


1.494868058150 


32 


0.095523702475 


14 


0.968484734380 


33 


0.000000000000 


15 


-0.945563574521 


34 


-0.879782299796 


16 


-1.101175502961 


35 


-1.833603623066 


17 


0.000000000000 


36 


0.234869257223 


18 


0.824250041644 
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Table 6: Fourier coefficients for / G M (ro(4), i, i?) 



R = 4.461438243496 R = 6.046497437542 



71 




dijl) 


Tl 




U 




f\ f\C\C\C\C\C\C\C\C\C\C\C\C\C\ 

U.UUUUUUUUUUUUUU 


U 


U.UUUUUUUUUUUU 


i 




1 nnnnnnnnnnnnnn 
i.UUUUUUUUUUUUUU 


1 


n nnnnnnnnnnnn 
U.UUUUUUUUUUUU 


Z 




U.Oooo4yoo44yUoD 


o 
z 


i .UUUUUUUUUUUU 


Q 

o 




n r: 1 7QQon/i 7/1 no 
U.OoOi ( oozy4 ( 4Uz 


6 


1 77070 rLQ^iQQ 71 

i. ( ( U ( yooDoo I i 


A 

4 




-U. ( U < iUD / oiioDD / 


A 

4 


u.uuuuuuuuuuzy 


r 





-U.UUUUUUUUUUUUUo 


r 




1 QQ 1 /1 70/1 0/i nriQ 
i.ooi4 ('U4y4UUo 







i.ZoyjoD (^UD4UUi4z 


D 


O 7/inQniXc:07i^7/1 

— U. / 4yoyooU r D /4 






u.yu ^ oDzooyiDDyo 


( 


n (17/1 Qf^OQ 1 Q^^OI^ 

U.U ^ 4oDyoioDzO 


Q 

o 




-U.44 ^ o4oooD moot) 


Q 

o 


n 70710(^701000 

U. / U / iUD ^ oizUy 


y 




n I^Of^ /I QQ70f^ /I Q77(^ 


y 


O 00000000001 

— U.UUUUUUUUUUlo 


iU 




n c;77(^Q/iOQn(^^?n70 

u.o/ mo4yoyDDy ^ z 


iU 


1 01/1/1i^1Q(^7QQ0 

— i.zi44oioD ( Ooz 


li 




1 1 0/1 Q IXQ77m IXi^ /I 1 

i.iz4ooo^ /yioD4i 


ii 


n i^007IXf^0/1 QQQQ 

U.DzU ^ oDz4oooo 


1 o 




-U.44yiooyU4Uor5oy 


1 o 
iz 


1 Oc:01 /1 1 7(^QOO/1 
i.z0zi4i / Dr5zU4 






U.UUUUUUUUUUUUUi 


io 


1 1OQ(^Q(^OO1i^0^? 
— i.iZrjDoDUziOOD 


1 A 




U.4oU (oy) ( looooZ ( 


1 /I 
i4 


O 71 10Q(^0Q1(^Q0 
— U. ( iizoDUoiOrjU 


io 




1 c:Qc;Qnmonnc;7Q/i 
-i.ooOrjyUizUUo / o4 


io 


O 0(^/1 70(^0Q0/1 1 Q 

U.yD4 / UDUoz4io 


Id 




U.oUUUUUUUUUUUio 


io 


O 0000000001 IO 

— U.UUUUUUUUUiiz 


1 { 




n 1 /I QQorii^noi /I /I QQ 
-U.i4oozUDyzi44or5 


1 7 
i / 


O 00000000010Q 

— U.UUUUUUUUUiUo 


1 C 




1 n(^/i7/innooni^(^ni^ 
i.UD4 ^ 4yUzzyoDUo 


io 


O 1 OQ(^/1 /I /I Q /I ^?00 

— U.izoD444o4DUy 


1 o 

iy 




U.zizDoy IDoDoDoZ 


iy 


/100(^(^/1Q71 1^0/1 

— u.4yyDD4o ( ioy4 


20 




0.00000000000003 


20 


0.941491814843 


21 




0.00000000000002 


21 


-0.447011950327 


22 




-1.56248056455209 


22 


-0.805501416928 


23 




0.85478501960318 


23 


1.236716356721 


24 




-0.90534916229569 


24 


-0.529902643583 


25 




0.45696099733973 


25 


0.000000008756 


36 




-0.37224839333969 


36 


-0.000000000168 


|c(4)c(9) 


- c(36)| = 4£;- 12 


288 


-0.064322242377 



|c(21)| = 2£; - 14 |c(2 • 32)c(2 • 4^) - c(2 • 122)| = 8£ - 11 

|c(17)| = 1£;- 10 
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Table 7: Supplemental table of Fourier coefficients for M (ro(2), 0, R) 



R 


8.92287648699174 


R 


12.0929948750786 


n 


A{n) 


n 


A{n) 


2 


-0 70710678118654 


2 


70710678118655 


3 


1.10378899562734 


3 


-0.70599475399569 


4 


0.49999999999993 


4 


0.49999999999999 


5 


0.90417459283958 


5 


-0.79974825694039 


6 


-0.78049668380711 


6 


-0.49921367803249 


7 


0.82934246755499 


7 


-1.71337067862845 


o 

8 


A O C O C C: O O AA C AO O A 

—0.35355339059330 


8 


A c c: c AA c: AOO 

0.35355339059328 


9 


0.21835014o8o776 


9 


A C:A1 CTI /IATOOACT/I 

—0.50157140733054 


10 


A r* O AO /I TAO C: ATO O A 

—0.63934798597339 


1 A 

10 


—0.56550741572467 


R 


13.77975135189073 


R 


13.77975135189073 




(even wrt z i— *■ — 77-) 




(odd wrt 2; 1-^ —77-) 


n 


A{n) 


n 


A(n) 


2 


2.96351804031448 


2 


13509091556820 


o 








4 


3.59139177031902 


4 


-0.79070303958101 


5 


0.73706038534834 


5 


0.73706038534830 


6 


0.73169192981688 


6 


0.03335391631439 


7 


-0.26142007576500 


7 


-0.26142007576538 


8 


2.60064131148226 


8 


-1.36013067551284 


9 


-0.93904050235826 


9 


-0.93904050236089 


10 


2.18429174878080 


10 


0.09957016228575 
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Table 8: Comparison of Fourier coefficients A{n) computed directly for M (ro(A^), 0, i?) {N 
1, 2), vs A{n) computed on M (ro(4), i, ^R) and using ijeTTI) . 



n 


A{n) 


Ain) 


\A{n) - 


A{n)\ 




feM{T 


o(2), 0,8.922876486992) {t = 


1) 




2 


"0.70710678118665 - 


-0.707106781186 


0.6£; - 


12 


3 


1.10378899562739 


1.103788995627 


0.4£; - 


12 


5 


0.90417459283969 


0.904174592840 


0.3£' - 


12 




/ C M(ro(l),0, 13.779751351891) (t = 


= 1) 




2 


1.54930447794126 


1.549304477941 


0.7^;- 


12 


3 


0.24689977245398 


0.246899772454 


0.3£;- 


12 


4 


1.40034436536892 


1.400344365369 


0.2£;- 


12 


5 


0.73706038534387 


0.737060385348 


0.4^;- 


11 


6 


0.38252292109716 


0.382522923066 


0.2£;- 


08 




/ € M(ro(l), 0, 13.779751351891) (t = 


:2) 




3 


0.24689977245437 


0.246899772454 


0.8£; - 


13 


5 


0.73706038535004 


0.737060385348 


0.2£;- 


11 


7 


-0.26142007624377 - 


-0.261420075765 


0.5£;- 


09 


9 


-0.93904050238904 - 


-0.939040502362 


0.3£;- 


10 



/ e M(ro(2), 0,12.092994875079) {t = 2) 
3 -0.70599475379863 -0.705994753996 0.2£; - 09 

5 -0.79974825694696 -0.799748256940 0.7i; - 11 

7 -1.71337067860377 -1.713370678628 0.2^; - 10 

9 -0.50157140750090 -0.501571407330 0.2^; - 09 

In segment 3, the calculation is based on (|6.ip and the second portion of Tabled! 
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Table 9: Comparison of Fourier coefficients for weigfits k = 9.044605824i? — 08 and fc = near 
an "avoided crossing". 



Corresponds to tlie cusp form 








k 


9.0446058240£;-08 





Difference 


R 


13.77975135189074 


13.77975135189074 






c(2) 


1.54930480559976 


1.54930447794069 


0.3^; 


-06 


c(3) 


0.24689988546553 


0.24689977245411 


O.IE 


-06 


c(4) 


1.40034433555250 


1.40034436536841 


O.IE 


-06 


c(5) 


0.73706067260516 


0.73706038534787 


0.2E 


-06 


c(6) 


0.38252272069428 


0.38252292306557 


0.2E 


-06 



Corresponds to the Eisenstein series 



k 


9.0446058240£:-08 





Difference 


R 


13.77975135138225 


13.77975135138225 






c(2) 


-2.06525760334129 


-1.98398933080188 


0.8E 


-01 


c(3) 


-1.72891679536648 


-1.68449330640991 


OAE 


-01 


c(4) 


2.97153986917404 


2.93621366473571 


OAE 


-01 


c(5) 


-2.02747287754385 


-1.96531634618530 


0.6E 


-01 


c(6) 


3.41008729668221 


3.34201674772446 


0.7E 


-01 



Tlie coefficients for the Eisenstein series at fc = were computed using (|5.ip . i.e.: 
c(2) = 2cos(i?ln2) 
c(3) = 2cos(i?ln3) 
c(4) = 1 + 2cos(i?ln4) 
c(5) = 2cos(i?ln5) 

c(6) = 2 cos(i? hi 6) + 2 cos(i?,(hi 3 - hi 2)) 

Table 10: Comparison of Fourier coefficients for weights k = 9.044605824E' — 08 and fc = "far" 
from an "avoided crossing". The weight coefficients were computed using the formulas in Table 

m 



k 


9.0446058240^-08 





Difference 


R 


13.62696884857618 


13.62696884857618 






c(2) 


-1.99957085683552 


-1.99957081810438 


0.4^;- 


07 


c(3) 


-1.48069687587703 


-1.48069680342062 


0.7^;- 


07 


c(4) 


2.99828354611637 


2.99828345661464 


0.7^;- 


07 


c(5) 


-1.99647405201235 


-1.99647406885962 


0.2^;- 


07 


c(6) 


2.96075820067617 


2.96075811858031 


0.8£:- 


07 



